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A  wind-tunnel  investigation  was  carried  out  on  an  oscillating  model  of  the  F-S 
winy  with  and  without  an  external  store  (A1M-9J  missile).  Detailed  steady  and 
unsteady  pressure  distributions  were  obtained  over  the  wing,  while  on  the  store 
total  aerodynamic  loads  were  measured. 

The  tests  covered  the  Mach  number  range  between  Ma  *  0.6  and  Ma  -  1.3S,  and 
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FOREWORD 


This  report  was  prepared  by  the  National  Aerospace  Laboratory 
(NLR  of  Amsterdam,  the  Netherlands.  The  sponsors  were  the  Air  Force 
Armament  Test  Laboratory  (A^ATL/DLJ)  of  Eglin  Air  Force  Base,  Florida 
and  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL/FBR  and  AFFDL/FBE) 
of  Wright-Patterson  Air  Force  Base,  Ohio.  The  sponsorship  was  per¬ 
formed  through  AFOSR  Grant  77-3233,  administered  by  Captain  D.  Wilkins 
of  the  Air  Force  Office  of  Scientific  Research  (AFOSR/TKN)  of  Bolling 
Air  Force  Base,  Washington  D.C. 

The  report  consists  of  four  parts.  Part  1  contains  the  general 
description;  Part  II  discusses  the  steady  and  unsteady  aerodynamic 
tests  of  the  clean  F-5  wing;  Part  III  discusses  the  tests  for  the  wing 
with  tip-mounted  stores;  and  Part  IV  discusses  the  tests  for  the  wing 
with  under-wing  stores. 

The  principal  investigators  were  Dr.  H.  Tijdeman  and  Mr.  J.  W.  G 
van  Nunen  of  NLR.  They  were  assisted  by  A.  N.  Kraan,  A.  J.  Persoon, 

R.  Poestkoke,  Dr  R.  Roos,  P.  Schippers  and  C.  M.  Siebert  of  NLR. 

Within  the  United  States  Air  Force,  this  program  was  initiated  by 
Lovic  Thomas  of  the  AFATL.  It  would  not  have  been  possible  without  the 
expert  assistance  of  Richard  Wallace  (Lt  Colonel,  USAF,  Retired),  and 
Lt  Colonel  Daniel  Seger  and  Major  Robert  Powell  of  the  European  Office 
of  Aerospace  Research  Development  (EOARD). 
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INTRODUCTION 


In  order  to  determine  the  unsteady  airloads  on  a  representative 
fighter-type  wing  at  transonic  and  supersonic  speeds,  wind-tunnel  tests 
were  carried  out  on  an  oscillating  model  of  the  F-5  wing  with  and  with¬ 
out  external  store.  If  present,  this  external  store  (AIM-9J  missile  + 
launcher)  was  mounted  either  at  the  tip  or  at  a  pylon  under  the  wing. 

The  wing  model  was  oscillated  in  pitch  about  a  50  per  cent  root- 
chord  axis  at  frequencies  varying  up  to  Uo  Hz  (for  dimensions  see 
Fig.  la).  The  Mach  number  ranged  between  0.6  and  1.35.  Detailed  pressure 
distributions,  both  steady  and  unsteady, were  measured  over  the  wing, 
while  on  the  store  the  total  aerodynamic  loads  were  obtained.  A  des¬ 
cription  of  the  experimental  test  set-up  and  the  test  program  is  given 
in  part  I  of  this  report  (Ref.  l).  The  results  are  published  in  a  data 
report  (Ref.  2),  while  for  easy  data  handling  they  are  available  also  on 
magnetic  tape. 

To  assist  in  the  evaluation  of  the  data,  reference  1  is  supplemented 
by  three  additional  parts,  covering  successively  the  clean  wing,  the  wing 
with  tipstore  and  the  wing  with  underving-store.  Each  part  contains  plots 
of  the  steady  and  unsteady  pressure  distributions  and  gives  a  brief  ana¬ 
lysis  of  some  selected  results. 

The  present  report  ( part  IV)  focusses  on  the  configuration  of  the 
wing  with  the  underwing-store  (figure  1b  shows  the  model  of  this  confi¬ 
guration  mounted  on  the  side  wall  of  the  NLR  High  Speed  Tunnel  (HST), 
while  the  location  of  the  pressure  orifices  is  given  in  figure  1c).  The 
report  considers  the  influence  of  the  pylon  and  store  on  the  wing  loading 
and  on  the  loads  experienced  by  the  store  itself.  For  the  wing  this  ana¬ 
lysis  concerns  a  comparison  with  the  clean  wing  configuration,  results  of 
which  are  taken  from  part  II  (Ref.  3).  For  the  loads  on  the  pylon  and 
store  the  contributions  due  to  the  pylon,  the  launcher,  and  the  missile 
body  with  aft  wings  and  canard  fins  are  considered.  This  last  part  of 
the  analysis  is  slightly  hampered  of  the  fact  that  during  one  of  the 
testruns  on  the  complete  configuration  the  canard  fins  broke  from  the 
model,  causing  the  omission  of  supersonic  data  for  the  complete  confi¬ 
guration. 

For  Ma  =  0.6  the  experimental  data  are  supplemented  with  theoreti¬ 
cal  results  obtained  with  the  NLRI -method,  an  unsteady  panel  method  for 
wing-body  configurations  (Ref.  5).  Further,  since  the  NLRI-method  is 


rather  -os.iy  in  terms  of  computer  costs,  a  successful  attempt  has  been 
mate  to  devise  a  panel  distribution,  which  will  lead  to  comparative  re¬ 
sults  when  applied  in  a  Doublet  Lattice  ca 


STFADY  RESULTS 


General 

.  iring  the  wind-tunnel  experiment  steady  pressures  were  measured 
r  the  v  ng  at  incidences  of  -0.5*  0  and  0.5  degrees.  A  listing  of 
the  test  variables  as  well  as  a  complete  set  of  plots  of  the  resulting 
steady  pressure  distributions  in  the  eight  measuring  sections  is  gath¬ 
ered  in  Appendix  IV. B.  In  addition  to  the  pressures  on  the  wing  alec 
*:he  normal  load,  the  side  load,  the  pitching  moment  and  the  yawing 
moment,  experienced  by  the  store,  as  well  as  the  side  load  and  yawing 
moment  felt  by  the  pylon,  were  measured.  These  loads  were  measured  by 
two  independent  s'rain  gage  balances  located  at  the  interfaces  of  wing 
and  pylon  as  well  as  pylon  and  store.  The  exact  locations  are  given  in 
figure  la.  The  results  are  given  in  Appendix  IV. D. 

For  the  limiting  case  of  zero  frequency  so-called  "quasi-steady" 
results  car  be  obtained  by  considering  the  steady  results  for  small 
incidences  around  a  certain  mid-position.  In  the  present  invest igati.r 
"quasi-steady"  loads  for  a  mid-position  of  zero  degrees  were  obtained 
from  the  steady  data  taken  at  incidences  of  -0.5  and  0.5  degrees.  The 
corresponding  definitions  are  given  in  Appendix  IV. A. 

The  steady  and  quasi-steady  load  distributions  on  the  wing  were 
obtained  form  the  measured  pressure  distributions.  In  this  respect,  it 
is  noted  that  for  the  integration  in  section  3  and  5  the  faulty  value, 
found  on  the  upperside  near  the  reading  edge,  was  replaced  by  a  new 
value  obtained  by  spanwise  interpolation  between  the  sections  2,  1 
and  6. 

Steady  and  quasi -steady  loads  on  the  wing 

In  order  to  better  understand  the  influence  the  undervir.g- store 
has  on  the  quasi-steady  and  unsteady  span vise  load  sidtributici  s ,  firs: 
its  effect  on  the  steady  load  distribution  will  be  considered.  Figures 
2a  through  >c  show  these  distributions  for  subsonic  (Ma  *  O.o',  tr&nscr- 
(Ma  *  O.y)  a..d  supersonic  (Ma  1.35)  flow  conditions.  Fo:  incac  ncec 
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of  -0.5,  0  and  0.5  degrees,  each  figure  gives  a  comparison  of  the  steady 
load  distributions  for  the  configuration  with  the  undervir.g-store  and 
for  the  clean  wing. 

In  general,  attaching  a  pylon  and  store  under  a  wing  has  a  consider¬ 
able  effect  on  the  flow  field.  The  flow  has  to  diverge,  which  locally 
results  in  increased  flow  velocities,  causing  a  further  lowering  of  the 
pressure,  especially  on  the  lower  side  of  the  wing  (e.g.  see  Figs. 

IV. B.  33  -  35  of  Appendix  IV.B.). 

In  figure  2a  it  is  shown  how  this  effect  influences  the  spanwise 
normal  load  distribution  at  Ma  =  0.6.  Besides  a  decrease  of  the  load, 
which  is  maximal  just  inboard  of  the  pylon  one  observes  also  a  discon¬ 
tinuity  across  the  pylon.  This  jump  in  the  loading,  which  corresponds 
to  a  side  load  on  the  pylon,  reduces  with  increasing  incidence. 

At  transonic  conditions  (Ma  =  0.9;  Fig.  2b)  the  addition  of  the 
pylon  and  store  introduces  a  region  of  supersonic  flow  over  the  lower 
surface  of  the  wing  on  both  sides  of  the  pylon.  This  supersonic  region 
is  responsible  for  the  relatively  large  decrease  of  the  steady  normal 
load  on  the  inboard  part  of  the  wing.  With  the  load  on  the  outboard  section 
remaining  almost  unaffected  by  the  presence  of  the  pylon/store  combina¬ 
tion,  the  jump  over  the  pylon  has  increased  with  respect  to  the  Ma  =  0.6 
case.  The  dependence  or  incidence  is  similar  to  the  subsonic  case,  be  it 
that  for  a  =  -0.5  the  normal  force  near  the  tip  reaches  even  a  larger 
value  than  found  for  the  clean  wing  configuration. 

For  Ma  =  1.35  (Fig.  2c)  the  interference  effect  from  the  pylon  and 
store  is  smaller  and  also  less  regular  over  the  span.  Just  inboard  of 
the  pylon  the  normal  load  is  decreased,  but  around  50  per  cent  of  the 
span  the  load  for  the  configuration  with  store  grows  above  the  value 
for  the  clean  wing.  In  section  2,  just  outside  the  zone  of  influence 
originating  from  the  front  of  the  pylon,  the  difference  becomes  negli¬ 
gible  again.  Between  the  pylon  and  the  wing  tip  the  interference  effects 
vary  strongly  with  the  spanwise  position.  The  jump  in  the  normal  force 
over  the  pylon  remains  small.  The  trend  with  incidence  is  similar  to  the 
other  cases. 

The  quasi-steady  spanwise  load  distributions  derived  from  the  steady 
distributions,  are  presented  in  figures  3a  to  3c.  Qualitatively,  the 
effect  of  the  pylon  and  store  is  the  same  for  the  three  flow  conditions 
considered:  an  increase  of  the  quasi-steady  loading  inboard  of  the  pylon 
and  a  decrease  on  the  outboard  part  of  the  wing,  together  resulting  in  a 
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tap  ol  iy  normal  force  across  the  pylon.  This  jump  is 

maximal  for  transonic  conditions,  for  which  the  pylon  and  store  have  i 

introduced  •.  r a  of  supersoni  •  flow  and  thus  a  shockwave  over  the 
lower  .surface  f  the  wing  on  bi  th  sides  of  the  pylon. 

To  •  f  feet  of  Mac:,  nu er  on  the  sectional  loads,  f '  gurt-  -t 

wive  th-  toasi -steady  normal  . oads  in  the  sections  2,  5  and  5  as  a 

wever,  .•  t  th<  parcitj  f  •  -  •  ■ 

for  th-3  ■  r.ri  durations  with  ur.derwing- store  a  clear  development  is  not 
■’  tn  I.  it  ■  :in  -  deluded  only  that  far  away  from  the  pylon  sections 
i::d  8)  thi  i::t*»rferene«  f  feet  s  are  negligible  for  all  Mach  numbers 

considered. 

.  toady  and  unsi-steady  Loads  i  the  pylon  ana  store 

preset  t<  t  -  ......  .  .  .  mderving-s  r«  Li 

various  stages  cf  completeness,  making  it  possible  to  investigate  the 
aerodynamic  contribution  of  the  pylon  and  of  specific  parts  of  store. 

The  behaviour  of  the  individual  forces  and  moments  with  Mach  number  is 
given  in  Appendix  IV. D.  In  figures  5a  and  5b  the  total  steady  side  force 
and  yawing  moment  experienced  by  the  combination  of  pylon  anc  store,  if 
present,  are  pi  -tted  versus  Mach  number  for  the  three  test  incidences 
-0.5  ,  0  and  0.5  degrees. 

v-  all  configurations  the  side  force  (Fig.  5a)  points  inward , 
while  the  magnitude  decreases  with  increasing  incidence.  Of  course, 
both  observations  are  in  accordance  with  the  behaviour  of  the  jump  in 
the  steady  spanwise  normal  load  distributions  dest ribed  in  sec. ion  2.2. 

Further,  for  all  incidences  the  addition  of  the  launcher  doubles  the 
total  side  force.  The  wit]  aft  win gs  and  the  canard 

has  a  simi  lar  effect,  although  due  to  the  mishap  with  the  canard  fins  no 
differentiation  can  be  made  for  the  effect  of  either  one  of  these  control 
surfaces . 

The  total  steady  yawing  moment  (Fig.  5b)  is  relatively  small  for 
ill  configurations,  and  no  clear  variations  with  incidence  are  apparent. 

The  contribution  of  the  pylon  i-  virtually  nil  for  all  cases  ‘us-,  ui 
^se-'  als  Fig.  IV. D. 6\  For  subsonic  Mach  numbers  the  launcher  does  not 
contribute  to  the  yawing  moment,  while  for  supersonic  conditio:."  *  he 
contribution  remains  small.  Most  of  the  yawing  moment  is  produced  by  the 
aft  wings  and  the  canard  fins. 


The  quasi-steady  side  force  and  yawing  moment  coefficients  for  the 
combination  of  the  pylon  and  store  are  presented  in  figure  6.  The  quasi¬ 
steady  side  force  is  positive  (pointing  outward)  for  all  configurations, 
which  is  in  line  with  the  behaviour  of  the  quasi-steady  spanwise  load 
distributions  on  the  wing  near  the  pylon  (Fig.  3).  With  only  the  pylon 
present, the  side  load  is  very  small.  Addition  of  the  launcher  increases 
the  side  load  uniformly  over  the  Mach  number  range  considered.  Next, 
adding  the  missile  body  together  with  the  aft  wings  and  canard  fins  more 
than  doubles  the  quasi-steady  side  load.  The  quasi-steady  yawing  moment 
is  small  for  all  configurations. 


3  UNSTEADY  RESULTS 

3.1  Vibration  modes 

To  monitor  the  vibration  modes  of  the  wing  and  store  during  the 
tests,  the  model  was  equipped  with  12  accelerometers,  eight  in  the  wing 
and  four  in  the  launcher,  of  which  two  were  measuring  in  vertical  direc¬ 
tion  and  the  other  two  in  lateral  direction  (for  exact  loacations  see 
figure  la).  In  the  same  way  as  for  the  clean  wing  configurations  (see 
Ref.  3)  also  for  the  configurations  with  the  underwing-store,  the 
readings  of  these  accelerometers  were  used  to  make  analytical  approxima¬ 
tions  of  the  vibration  modes  for  the  20  Hz  testruns.  The  polynomial 
expression  used  for  these  approximations  is: 

/  ,  v  n  m  (m=0,  1,2;  n=0,1  for  the  wing) 

w(x,y)  =  l  a^x  y 

(m  =  0;  n  =  0,1  for  the  store) 

This  expression  assumes  no  deformation  in  chordwise  direction  and  a 
parabolic  deformation  in  spanwise  direction.  In  using  the  measured 
values,  for  testrun  1*8  the  reading  of  accelerometer  5  was  weighted  to 
ensure  a  realistic  development  of  the  torsion  angle  along  the  span. 
Further,  for  testrun  1*3  the  accelerometers  in  the  store  measuring  in  ver¬ 
tical  direction  gave  erroneous  readings.  Therefore,  the  vertical  dis¬ 
placement  of  the  store  was  taken  equal  to  the  displacement  of  the  wing 
at  the  location  of  the  pylon  attachment ,  which  was  in  accordance  with 
the  findings  for  the  other  testruns.  Taking  these  two  points  into  account, 
the  difference  between  the  analytically  approximated  displacements  and 
the  experimental  readings  remains  within  2  per  cent  of  the  maximal  dis¬ 
placement.  Table  1  presents  the  numerical  values  of  the  coefficients  a 
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for  *ln  1:1m.  \.im  vibration  modes .  The  normalization  is  carried  out  such 
that  at  tin'  v  in*  section  of  accelerometer  the  tangent  of  the  angle  of 
oscillation  equals  one. 

To  obtain  insight  into  the  effect  of  the  pylon  and  store  c*s  the  v;b ra¬ 
tion  modes  of  the  i  n  k  .  the  nodal  linen  us  found  for  the  clean  w  ng  (see 
ri  fence  t  iuid  the  wing  with  underv .  ug-st.ore  are  presented  n  ;  gures 
"a  a:  i  h.  for  Ma  J.fc  and  i'.  ■>  the  comparison  concerns  the  store  with 
the  canard  ‘‘ins  (Fig.  'a',  while  for  Ma  *  0.9,  1.10  and  1 .  tr'  the  store 
t 'iout  the  canard  f;ns  v  h'  i  l' '  is  considered,  ''.'tie  figures  show  that 
for  .'0  II.-  1  c  addition  of  the  pylon  store  hardly  influences  the  tvs  it.  ion 
of  the  nods  1  ie.  The  difference  found  for  Ma  ■  1 .  1 0  (fig.  ’b)  most 
probably  i  caused  y  ‘he  change  in  stagnation  pressure  from 
'  * ' 1  v  '  '  for  the  clean  wing  ‘o  0.7  x  'O'  Pa  for  the  wing  with  under¬ 
wing-store.  part  111  of  this  report  (Kef.  1« '  it  was  found  that  for  the 
wing  with  t  pat  ore  such  a  change  in  stagnation  pressure  makes  the  nodal 
:nc  move  lVrvarl,  which  corresponds  with  the  behaviour  for  the  present 
configuration . 

The  development  with  Mach  number  is  similar  to  that  of  the  clean 
wing  i.o.:  in  the  subsonic  and  transonic  regime  the  nodal  line  bends 
to  the  rear,  while  in  the  supersonic  range  the  nodal  line  tends  to  re¬ 
turn  to  its  original  position.  The  local  torsion  angle  was  found  to 
remain  almost,  constant  over  the  wing.  As  mentioned  before,  in  vertical 
direction  the  store  follows  the  motion  of  the  wing.  In  lateral  direction 
the  motion  is  negligibly  small. 

No  appj  <xi  mat.  ions  wen-  made  for  the  V  He.  test  runs  of  the  configu¬ 
ration  with  nnderw ing- rt  ore.  However,  s  preliminary  analysis  wa  msec 
check  lie  validity  of  the  results.  Contrary  to  what  was  found  for 
the  wing  with  complete  tips  to  re  (Ref.  •'* )  all  runs  except  no.  show 
•i  stable  vibrat  ion  mode.  In  test  run  a  complex  mode  was  euccutit ,  vd . 
which  casts  doubts  on  the  unsteady  store  loads  measured  by  the  balancer. . 


.'Udy  l  "Ids 

Journal 

1  r  all  pylon/ store  configurations  the  unsteady  pressure  1  >  men¬ 
tions  measured  over  the  wtng  are  tabulated  in  reference  and  j  t ;  ed  n 
the  present  -aport.  I'he s  •  plots  and  a  listing  of  the  releva-r.  i 
-» r  i  -a  bit  -Vi-  -  e  m  A,  pcjulix  J.C.  For  a  better  compart  on 
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values 


on  the  upperside  are  shown  with  u  reversed  sign. 

In  the  present  experiment  the  unsteady  side  force  and  yawing 
moment  acting  on  the  pylon, as  well  as  the  unsteady  normal  force,  side 
force,  pitching  moment  and  yawing  moment  on  the  store  were  measured  by 
the  same  strain  gage  balances  used  to  obtain  the  steady  data  (for  details 
see  references  1  and  6).  Plots  of  these  quantities  versus  Mach  number 
are  presented  in  Appendix  IV. D  for  all  20  Hz  test  oases. 

The  problems  with  the  vibration  modes  in  the  test  runs  43  and  U U 
are  the  cause  for  the  incorrectness  of  the  values  for  the  forces  find 
moments  on  pylon  and  store  as  given  in  the  tabulated  results  (Kef.  2). 

For  run  U 3  the  correct  values  could  be  retrieved.  They  are  presented  in 
table  2. 

3.2.2  Unsteady  loads  on  the  wing 

The  development  of  the  unsteady  normal  force*  with  Mach  number  in 
the  sections  2,  5  and  8  is  given  in  figure  8.  For  a  frequency  of  20  Hz 
this  figure  presents  a  comparison  between  the  results  of  the  clean  wing 
and  the  wing  with  underwing-store.  For  the  lower  Mach  numbers  the  store 

is  complete,  while  for  the  higher  ones  the  data  represent  the  configura¬ 

tion  with  the  canard  fins  missing.  At  Ma  «  0.9  both  configurations  were 
tested. 

Comparison  of  figures  4  and  8  show  that  the  character  of  the  inter¬ 
ference  found  for  20  Hz  is  similar  to  that  for  quasi-steady  conditions 
except  in  section  8.  In  section  2  the  presence  of  the  pylon  and  store 
is  not  felt,  while  in  section  5  the  difference  between  the  two  configu¬ 
rations  at  20  Hz  is  slightly  less  than  for  0  Hz.  However,  in  section  8, 

contrary  to  the  quasi-steady  findings,  the  pylon  and  store  do  influence 

the  unsteady  normal  load  at  20  Hz.  As  explained  below,  the  reason  for 
this  is  different  for  each  Mach  number.  The  large  differences  found  for 
Ma  =  1.1  may  be  caused  partly  by  the  change  in  stagnation  pressure  from 

P  =  1.0  x  105  Pa  to  P  =  0.7  x  105  Pa. 
o  o 

To  find  out  how  the  interference  of  the  pylon  and  store  spreads  over 
the  wing,  in  figures  9  to  1 1  the  spanwise  distribution  of  the  unsteady 
normal  force  and  pitching  moment  for  the  configuration  with  and  without 

*) 

For  the  integration  of  the  unsteady  pressure  distributions  in  sections 
3  and  (j  the  zero  value,  wrongly  measured  on  the  upperside  near  the 
leading  edge,  was  replaced  by  a  new  value,  obtained  by  spanwise  inter¬ 
polation  between  sections  2,  4  and  6. 
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store  are  conparcd  for  Ma  *  0.6,  0.0  and  t .  3S  mi.,!  «  freouer  oy  f 

At.  Ma  »  0 . (>  t!ie  pylon  and  store  introduce  the  expert*  '  •••  ;•  the 

real  part  .  f  the  normal  force  distribution  ( Fig.  ■>«) .  however,  contrary 
to  tile  quasi-stead}  case  (Fig.  \a ) ,  this  interference  retaa '  •  >-d 

around  the  attachment  position  of  the  pylon.  The  inboard  V  per  vis 
spar  apparent Ly  is  not  affected  by  the  rresence  of  the  store. 
reason  found  by  comparing  the  unsteady  pressure  tistribut  inn  for  the 
two  configurations.  This  shows  that  for  the  wing  with  st  or--  • ' 
tend}-  pressures  at  the  Leading  edge  on  the  lower  side  are  s  •' gu f  leant 
lower  than  for  t  he  clean  wing,  thus  resulting  in  a  lover  unst  e.-t  iy  norm 
force.  The  relat  ively  large  difference  found  in  section  Is  oar.  be  -  x- 
p  laired  by  t  lie  fact  that  the  bulge  found  in  the  unsteady  pressure  dist-r 
but  i  -u  on  the  lower  side  of  the  clean  wing  (sop  part  T1  of  t  h.  present 
report,  Hof.  V,  does  not  show  up  for  the  present  configurn'  .on  (see 
V:  g.  IV.  0.  j.’  of  Appendix  IV.  C) .  The  imaginary  part  of  the  unsteady  nor¬ 
mal  force  is  not  influenced  at  all  hy  the  pylon  and  store.  As  'ar  ns  the 
pitching  moment  distribution  is  concerned  (Fig.  <5b) ,  the  real  part  shows 
an  interference  more  spread  out  over  the  wing;  however,  it  must  be  re¬ 
alised  that  the  values  remain  very  small.  The  difference  between  t'e  t wo 
conf Igurations  corresponds  to  a  shift  of  the  normal  force  of  a' out  ’■  per 
cent  of  the  local  chord.  As  with  the  normal  force,  the  imaginary  part  is 
free  of  interference  effects. 

At  transonic  flow  conditions  (Ma  *  ('.o'  the  pylon/store  interference 
extends  over  the  whole  wing.  Inboard  of  the  pylon  the  normal  force  dis¬ 
tribution  (Fig.  10a)  shows  an  irregular  interference  pattern.  The  shape 
*  '■  »  rer.l  nart.  of  tb  •  .-t  -  ndy  spa". wise  list  ri  but  ion  for  the  v  :u  v  5 1  *• 

.’.ore  -or responds  very  well  with  the  quasi-steady  one  (see  Fig.  (  the 
tip  at  section  }  is  caused  by  the  shock  which  for  this  conf igurat ion 
exists  or.  the  lower  wing  surface.  The  sharp  decrease  at  section  b  was 
exo  lined  in  part  II  (Kef.  1'  as  partly  due  to  an  insufficient  e-'.-liit  on 

■  mry  1  :•••  .  jr.-it  '  over  tin  peaks,  generat  ed  by  the  h--e  -\  •••• 

-  u'p  r  .  fa  ■  the  wing.  Outboard  of  the  pylon,  the  presen  >e  -f  Ms- 

1  gni  fi  cant  deersase .  1  .  -i-  lary  part  f  tl  icr 

mal  load  i  •  af fe -ted  most  near  the  wing-tip.  As  can  be  exp-'ct*  1  f' 
t-.t  -hmg  m  flue-iced,  m  '.  so  rather  irregularly  ivig.  HT-  .  i'he 

maximum  in'  e  -ference  is  fo  unt  in  the  vicinity  of  the  pylon. 

l.i  only  neasu’v- !•>:.'  for  Mil  '  i\'>  the  '  >nf  gun- 

u:ido-wai(  it.,  -re  was  test  e«t  with  t  without  canar  l  fins  in  or.hu 
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on  overlap  between  the  subsonic  and  supersonic  case,  of  which  the  latter 
was  tested  with  only  the  aft  wings  present.  The  figures  10a  and  10b  show 
that  the  influence  of  the  canard  fins  on  the  unsteady  spanwise  load  dis¬ 
tributions  is  negligible. 

At  Ma  »  1.35  the  effect  of  the  underwing-store  is  very  regular.  The 
real  part  of  the  normal  force  (Fig.  11a)  is  increased  over  the  inboard 
part  of  the  wing  and  decreased  over  the  outboard  part,  as  was  observed 
also  for  the  quasi-steady  case  (Fig.  3c).  In  section  2, just  outside  the 
zone  of  influence  originating  from  the  front  of  the  pylon,  the  influence 
has  become  almost  zero.  The  imaginary  part,  being  smaller  in  magnitude, 
is  affected  only  near  the  pylon.  The  behaviour  of  the  unsteady  pitching 
moment  is  similar  apart  from  a  dip  in  section  1*  (Fig.  1 1  b ) . 

3.2.3  Unsteady  loads  on  the  pylon  and  store 

The  sum  of  the  unsteady  lateral  loads  at  20  Hz  measured  by  the 
balances  in  the  pylon  and  the  store  is  plotted  in  figure  12  versus 
Mach  number.  As  the  configuration  is  gradually  made  more  complex,  the 
influence  of  the  different  parts  of  the  store  can  be  studied. 

The  behaviour  of  the  real  part,  of  the  side  force  and  the  yawing 
moment  coefficients  is  very  similar  to  what  was  found  for  0  Hz  (Fig.  6). 
This  is  valid  for  the  development  with  Mach  number  as  well  as  for  the 
variations  resulting  from  changes  in  the  configuration.  The  main  contri¬ 
bution  to  the  side  force  comes  from  the  aft.  wings.  For  Ma  =  0.9  it  is 

shown  that  the  canard  fins  give  a  negligible  contribution.  The  reason 
for  this  is  clear.  Their  location  far  in  front  of  the  main  wing  makes 
any  interference  originating  from  this  wing  minimal.  Therefore  the  only 
way  the  canard  fins  can  contribute  to  the  side  force  is  by  means  of  a 

lateral  motion  of  the  store  itself,  and  this  motion  was  found  to  be 

negligible  (see  section  3.1).  The  imaginary  part  of  the  side  force  is 
almost  zero,  with  only  a  small  contribution  of  the  aft  wings  in  the 
supersonic  range. 

The  total  unsteady  yawing  moment  remains  small.  As  fas  as  the  real 
part  is  concerned,  the  launcher  and  missile  body  contribute  most,  while 
the  aft  wings  give  a  small  reduction.  The  imaginary  part  is  negligible. 
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T r«  NLRI-ra-thod 

Next  to  the  Doublet  Lattice  method  the  NLP  possesses  a  theoretical 
net h  ;>l  for  the  calculation  of  ‘  msfrady  aerodynamic  forces  on  a 

»  e  flow.  This  so-called 

"NLhl-methi.  1”  ••  h  panel-type  ho.*  based  on  the  potential  flow  apprex- 

imat  ion.  F-  r  *  he  >dies,  such,  as  fuselage  and  stores,  use  is  made  of  a 
repr  •  ••n‘ a*  ion  •  •  nos  of  unsteady  source  panels  distributed  ?v-r  the 
contour  of  the  body,  while  for  the  lifting  surfaces,  which  are  taken 
to  he  infinitely  thin,  the  Doublet  Lattice  approach  is  used.  A  detailed 
description  of  the  method  .s  given  in  reference  % 

As  par*  f  the  present  investigation  the  NLRI-method  was  applied 

•  c  f  gural  .  with  the  coaplete  underwing-store  at  Ma  =  o.o  and 

a  frequency  of  20  Hz.  I’he  piuie-  distribution  used  in  this  calculation  in 
shown  in  figure  13.  The  wing,  the  pylon,  the  canard  fins  and  the  aft 
wings  are  represented  as  thin  lifting  surfaces,  while  *  he  launcher  and 
missile  body  are  combined  to  one  long  body  with  an  open  rear  end.  Apart 
from  the  pylon,  the  panel  distribution  i‘.  similar  to  the  one  us".  • 
the  NLRI-cal  dilations  made  for  the  wing  with  tiprtore  (see  l’art  Ml, 

Ref.  U ) . 

The  the'  retical  pressure  distributions  in  three  sect  Ions  of  the 

* )  . 

wing  arc  compared  with  experimental  data  in  figure  1u.  The  theoreccal 
lines  represent  the  calculated  pressure  jump  across  the  wing  divided  by 
‘  wi  .  In  he  sections  2  and  8  the  agreement,  between  theory  anu  experiment 

•  *  good.  In  section  ‘  ,  j\:  t  aboard  of  the  pylon  the  agreemetr  oks  a  ' 

mood.  However,  it  should  be  noted  that  for  this  section  because  f  the 
. ocai  asymmetry  cf  the  configuration  the  comparison  would  be  bet'  r  f 
the  experimental  data  would  represent  also  the  jump  across  the  wing  sur¬ 
face.  In  fact,  the  theoretical  line  quit<  closely  averages  the  experi¬ 
mental  values.  The  remaining  differences  may  he  explained  by  '  . . 

ii"  the  '’low  stagnates  agains'  the  pylon  (sec  Fig.  1V.B.  3-  o!  Appi  — 
nix  IV. B),  an  effect  which  cannot  be  accounted  for  in  the  theory. 

A  corr.rarison  of  the  corresponding  theoretical  and  experimet  t  a  I  un- 
ndy  sj unwise  normal  loau  distr  butions  for  the  wing  vi*h  an  v  • h  \  ' 

Note  that  for  the  purpose  of  comparison  the  unsteady  pres  to-  - 

taped  •  .  upper  Lde  are  plotted  with  a  wj  Lgn. 


undcrw lag-store  is  given  in  figure  r..  The  theoretical  results  for  the 
wing  with  store  were  comouted  with  the  NLKI-method,  while  for  the  clean 
wing  the  Doublet  Lattice  method  was  used  (Kef.  i) . 

As  found  for  the  clean  wing,  also  for  the  wing  with  pylon  and  store, 
the  real  part  of  the  calculated  unsteady  normal  force,  as  normal,  is 
larger  than  the  experimental  value.  Since  the  differences  are  similar 
for  both  configurations ,  conclusions  concerning  interference  effects 
should  remain  valid. 

The  theoretical  calculations  predict  an  interference  due  to  the 
presence  of  the  pylon  and  store,  which  agrees  reasonably  with  the  ex¬ 
perimental  findings.  In  the  real  part  of  the  unsteady  loading  the  calcu¬ 
lated  and  measured  jump  over  the  pylon  are  of  the  same  order.  The  pre¬ 
dicted  reduction  of  the  load  on  the  outboard  part  of  the  wing  is  larger 
than  found  experimentally.  The  reason  for  this  is  the  relatively  low 
experimental  values  for  the  clean  wing,  which  as  explained  in  part  II 
(Kef.  3)  are  caused  by  unexpected  peaks  and  bulges  in  the  measured  un¬ 
steady  pressure  distributions  and  which  tend  to  lower  the  integrated 
values.  The  small  increase  of  the  normal  load  inboard  of  the  pylon  is 
not  observed  in  the  tests  for  20  Kz  for  reasons  explained  above  in 
section  3.2.2,  but  a  similar  increase  is  found  in  the  quasi-steady  case 
(see  Fig.  3a). 

In  table  3  the  experimental  normal  force,  pitching  moment,  side 
force  and  yawing  moment  acting  on  the  combination  of  pylon  and  store 
are  compared  with  the  values  obtained  with  the  NLRI-method.  For  the 
unsteady  normal  and  side  forces  the  agreement  is  good,  while  for  the 
pitching  and  yawing  moment  the  theoretical  values  are  off  by  a  factor 
of  two.  However,  in  evaluating  this,  it  should  be  realized  that  the 
normal  force  as  well  as  both  moments  are  very  small.  For  example,  the 
real  part  of  the  normal  force  is  only  T  per  cent  of  the  force  acting 
on  a  strip  of  the  wing  on  both  sides  of  the  pylon  having  the  same  width 
as  the  store.  For  the  side  force  this  value  is  about  30  per  cent. 

Schematization  in  terms  of  the  Doublet  Lattice  Method 

Compared  to  the  Doublet.  Lattice  (PL)  method,  the  NLKI-method  uses 
2  to  \  times  more  computer  time.  This  makes  the  computation  of  unsteady 
aerodynamic  force  coefficients  for  flutter  certification  very  expensive 
and  impractical  if  carried  out.  solely  with  the  NLKI-method.  Therefore, 
as  proposed  in  reference  5,  it  is  worthwhile  to  investigate  the 


puss  It  *lity  to  . urry  out  such  calculations  solely  with  a  cheaper  method, 
such  as  the  DL-method,  after  the  panel  distribution  has  been  checked 
against  renal*:,  f  a  Nil!  I -computation. 

In  the  present  investigation  it  is  tried  to  devise  such  a  panel 
distribution  for  the  DL-meth.'d ,  in  which  the  store  is  represented  by 
additional  thin  ifting  surface  part s .  "hose  additional  lifting  surfaces 
have  to  :'u*fil  tv  roles.  They  should  lead  to  a  good  prediction  of  the 
interference  or.  the  wing  and  at  the  same  time  give  reliable  values  for 
he  to:  -oads.  In  part  III  (Ref.  k)  a  panel  distribution  is  shewn 
which  satisfies  these  criteria  reasonably  well  for  "he  configuration 
with  complete  tipstore.  The  same  panelling  scheme,  with  an  additional 
horizontal  row  of  panels  to  model  normal  force  and  pitching  moment  on 
the  launcher  and  missile  body,  is  used  for  the  wing  with  underwir.g- 
store.  This  panel  '•  ution  is  given  in  figure  16. 

The  ui  steady  pressure  and  normal  load  distribut  ions  on  the  wing, 
obtained  with  this  DL  panel-layout  for  Ma  =0.6  and  F  =  20  Hz,  are 
identical  to  the  ones  found  with  the  NLRl-nethod.  The  unsteady  loads  or. 
pylon  and  store  are  given  in  table  3.  Although  the  panel-layout  was 
adopted  from  the  configuration  with  tipstore  and  no  special  tuning  for 
the  present  configuration  was  made,  the  agreement  with  experimental  store 
loads  is  very  good.  Even  the  moments  are  predicted  well. 

Clearly,  the  type  of  panel  schematization  presented  above  for  the 
Doublet  Lattice  method  may  be  of  use  in  practical  flutter  investigations. 

CONCLUSIONS 

A  wind-tunnel  investigation  was  carried  out  on  a  harmor. : cal ly  oscil¬ 
lating  model  of  the  F-5  wing  with  and  without  external  store.  As  part  of 
the  test  program  steady  and  unsteady  measurements  were  per forme,.  '  the 
configuration  with  a  pylon  and  store  mounted  under  the  wing.  •  - 

rem.or.tr  concerned  pressures  or.  the  wing  and  aerodynamic  1  ■  ,  .  * !  • 

pylon  and  store.  The  present  report  gives  a  brief  analysis  of  the  result 
for  this  configuration.  This  analysis  concerns: 

.  the  interference  of  the  tipstore  on  the  steady,  quasi  ->.!.y 
and  unsteady  vir.g  loadings,  indicating  that  pylon  ur.e  re  are 
responsible  for  an  increase  of  the  loading  inboar  1  ‘  p\  .'■■■ 

a i.J  a  leorease  in  the  outboard  side 


the  behaviour  vith  Mach  number,  which  shows  that  at  Ma  *  0.9 
the  wing  loading  is  dominated  by  the  formation  of  a  Bhock  on 
the  lower  aide  of  the  wing 

the  influence  of  the  complexity  of  the  store,  which  indicates 
that  the  side  force  grows  consistently  as  the  store  becomes 
more  complex,  while  the  yawing  moment  remains  very  small 
a  description  in  terms  of  polynomials  of  the  in-wind  vibration 
modes  for  20  Hz 

a  comparison  of  these  vibration  modes  with  those  of  the  clean 
wing,  showing  that  the  pylon  and  store  hardly  influence  the  posi¬ 
tion  of  the  nodal  line 

a  comparison  of  unsteady  experimental  data  for  Ma  »  0.6  with 
the  results  obtained  with  the  NLRI-method,  showing  a  satis¬ 
factory  agreement 

the  definition  of  a  panel  distribution  for  Doublet  Lattice  cal¬ 
culations  of  which  the  results  agree  well  with  the  experimented, 
data  and  which  may  be  of  use  in  practical  flutter  investigations. 
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wind  tunnel 


Fig.  3a  Quasl-ateady  spanwise  normal  load  dint ribut Iona  on  the 
clean  wing  and  the  wing  with  pylon  and  store  at  Ma  «*  0. 
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Fig.  3a  Steady  side  force  on  the  pylon  and  store  for  various 
configurations 


-0.0021 


-0.002L 


1.4  Mo 


a  =  +o.5° 


+  □ 

Pq  =  1.0  *  105  Po 

■f  cr  A 

P0  =  0.7  X  105  Po 

1.4  Mo 


1.4  Ma 


a  *  -0.5° ! 


FI*  Sb  Steady  ynwing  rcomont  on  the  pylon  and  store  for  vari,.  .-- 
conf  lgurntlons 


F  -  0  Hi 

CONFIGURATION 


FYLON  (P) 

P  ♦  LAUNCHER  (  L  ) 

P  ♦  L  +•  MISSILE  BODY  (MB) 

+  AFT  WINGS  (AW) 

P  +  L+  MB+AW  +  CANARD  FINS 


Fig.  6  Quasi-steady  side  force  and  yawing  moment  on  the  pylon 
and  store  for  various  configurations 
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Doublet  Lattice  method 
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APPENDIX  IV. A. 

Definitions  of  Steady  and  Unsteady  Aerodynamic  Quantities* 


THE  WING 


Steady 


’’ressore  coefficient  C  : 

P 

Cp  *  (Ploc  -  • 

Sectional  normal  force  : 


Z  =  C 


zQC  •  Cz  =  -  /  'Cp+  -  CJ  d(x/C) 


oectional  pitching  moment  about  quarter  chord  point  (positive 


nose  down)  : 


M  =  C^QC 2 ,  Cm  =  -  /  (cp+  -  Cp_)  (x/C  -  0.25)  d(x/C) 


Quasi -steady  at  zero  incidence  (w  =  0:  n  =  0°) 

o 

Pressure  coefficient  C  : 

pq 


C  =  AC  /Aa 

pq  p 


Cp  (ao  +  Aa1>  -  C  (ao  -  Aa  ) 


Aa  +  Aa 
1  2 


Sectional  normal  force  : 


za  =  ffQCCzo9elWt  .  C  =  -  AC  /Aa  =  1 

q  zq  zq  it  z  7r 


=  I  Cz  (a0  +  Aal}  -  Cz  (cto  -  Aa2) 


Aa-j  +  Aa2 


Sectional  pitching  moment  (positive  nose  down)  : 

Mq  -  |  QC2C  9eiut  ,  C  =  -  AC  /Aa  =  -  -m>°  *  ^  ^  ~ 

q  ^  nq  mq  it  m  tt  Tia T~Ka 

1  2 


uiic  definitions  for  the  unsteady  aerodynamic  quantities  are  accc;  • 
to  the  AGAP.D  manual  on  Aeroelasticity  vol.  VI  ^Hef.  6). 


rnamammmm 
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i  ■ 


hunt  en.tj 


Pressure  cor t'f  i  o  i fiit  r 


p> 


C  .  -  HeC  ♦  ilmO 

pi  pi  pi 


r./Q 

1 


Poet.  iona  l  normal  force  : 


■  nQCC  .(!•“'*  ,  ('  ■  lid’  ♦  i  I  mi' 

i  7.  i  .-i  .-i 


/  (0  .  -  C  .  )  d(x/C) 

n  pi-*-  01- 


1"- 


feet ional  pitching  moment  about  quarter  chorit  point  (positive  none 
itown  )  : 


..  n  . 

M.  ■  i  Ql‘  i  (' ■' 

i  .  m  i 


mi 


Ret'  .  +  i  lmC  .  ■  -  -  f  (f  -  f  .  )  tx/c'  -  d(x/C) 

mi  mi  ti  1  pi+  pi- 


THK  :'\  U'N  Affli  ;ViMRK 

St  i-.-nh- 

Normal  fore<  : 

-  C..QOS  . 

U 

Side  force  (  po.-i  i  I  i  Vc  outward'  : 

Y  -  CyQCS  . 

Pitching  moment  about  balance  centre  l posit  iv<-  noae  up) 

M  -  l’  tJOaS  . 

M 

Yawing  moment  about  balance  centre  vponitive  none  inward'  : 
N  -  . 


APPENDIX  IV.B 


Steady  Pressure  Distributions 


*) 


Wing  +  Pylon  (Conf.  10) 


Note  that  in  sections  3  and  b  the  first  point,  on  the  upperside  shoving 
a  zero  value  is  a  faulty  pressure  point,  which  should  not  be  consider¬ 
ed  in  any  evaluation. 
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CONF.  31  (WING  +  PYLON  ♦  LAUNCHER  +  MISSILEBODY  WITH  AFf  WINGS) 

,5, 


Fig.  IV. D. 3 


STEADY  NORMAL  FORCE 
ON  THE  STORE 


CONFIGURATION _ 

o  PYLON  (P)+ LAUNCHER  (L) 

A  P  +  L  + MISSILE  BODY  (MB) 

+  AFT  WINGS  (AW) 

□  P  +  L  +  MB  +  AW  +  CANARD  FINS 


CT=  +0.5° 


O'  =  0° 


1.0  1.4  Mo 


□ 

P0  =  1.0  x  105  Po 

O'  A 

P0  =  0.7  x  105  Po 

|<>  =  -0.5° 


0.6 


1.0 


1.4  Mo 


CONFIGURATION 


PYLON  (P)  +  LAUNCHER  (L) 

P  +  L  + MISSILE  BODY  (MB) 

+  AFT  WINGS  (AW) 

P+L  +  MBt  AW+  CANARD  FINS 


UNSTEADY  SIDE  FORCE 
ON  THE  STORF 


UNSTEADY  YAWING  MOMENT 
ON  THE  STORE 


□  A 

P0  =  1.0  *  105  Po 

cr  a 

P0  =  0.7  *  105  Po 

